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Abstract 

A complete quantum field theoretic study of charged and neutral par- 
ticle creation in a rapidly/adiabatically expanding Friedman- Robertson- 
Walker metric for an 0(4) scalar field theory with quartic interactions 
(admitting a phase transition) is given. Quantization is carried out by in- 
clusion of quantum fluctuations. We show that the quantized Hamiltonian 
admits an su(l,l) invariance. The squeezing transformation diagonalizes 
the Hamiltonian and shows that the dynamical states are squeezed states. 
Allowing for different forms of the expansion parameter, we show how the 
neutral and charged particle production rates change as the expansion is 
rapid or adiabatic. The effects of the expansion rate versus the symmetry 
restoration rate on the squeezing parameter is shown. 



1 Introduction 

In the cosmological context, early studies on the creation of particles in curved 
spacetimes have concentrated on neutral scalar fields [1] , [2] . While these studies 
led to an understanding of relationships between particle physics in the early uni- 
verse and cosmology, the inclusion of spontaneous symmetry breaking coupled 
with the HEP scenario for the early universe has led to the recent resurgence in 
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inflation being driven by the fluctuations in quantum fields undergoing symme- 
try breaking (preheating) [3, 4]. The basic methodologies used for such studies 
has largely been the Landau -Ginzburg model [6]. In this paper we advocate 
the use of a Hamiltonian generated through an 0(4) linear sigma model with 
quartic interactions for such studies. The Hamiltonian is a generic Hamiltonian 
applicable to many domains [5]. In this paper we show how it can be used to 
determine the production rates of charged vs. neutral particles in an expanding 
metric. We also show how squeezed states are the correct dynamical states for 
the fields in expanding spacetimes. We hope to show in a later communication, 
how these states may be used to study the "tachyonic" inflationary cosmological 
models of current interest. 



2 The Hamiltonian. 



The growth of long wavelength fluctuations during breaking of symmetries has 
been studied recently for the case of a single scalar field with quartic interactions 
[3]. While they have used Lattice simulations, we advocate a model Hamiltonian 
for quantum fluctuations of charged and neutral scalar fields which we shall 
derive here. Our starting point is the 0(4) linear sigma model with an action 
given by: 



S = j d?xdta{tf{\i> 2 2?(V*0 a " V2f) 



(1) 



with 



(2) 



where i — 1..4 are real scalar fields. The background spacetime is taken to 
be a FRW model: 



ds 2 = dt 2 - a(t) 2 dx 2 , 



(3) 



where a(t) is the expansion parameter. We now use a background field analysis 
to study the quantum effects. Assume $j has a background classical component 
fa which satisfies the classical equations of motion: 



5S_ 
5$, 



0. 



Treat the quantum field fa as fluctuation around classical solution: 

— >fa + fa 

since fa satisfies the classical equations of motion, 



1 



S = S[fa] + - 



5 2 S 



T l 6<!>i6<f>. 



*=4 



+ 



(4) 



(5) 



(6) 
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We shall deal with a quadratic fluctuation action given simply by: 

„ 1 f 5 2 S r 



(7) 



For the particular scalar field action given above, assuming all fields vanish at 
infinity (asymptotically flat metrics), 



5S 



, SV 



Imposing the classical equations of motion, we find that 

BV 

O^g^d^i) + a 3 m 2 fa + a 3 — = 



(8) 



(9) 



where = j^-\<j H - The equations of motion in this metric are 



A 
a 



1 

a 



-^V 2 ^ + m 2 </>, + ^- = 0. 



dV 

84>i 



(10) 



Since we are interested in the dynamics of the fluctuation field, we shall treat 
the fluctuation field in S 2 as a classical field and S2 itself as the classical action 
for its dynamics. We define a Lagrangian density for studying the dynamics of 
the fluctuations, £, as follows: 



2 2 -^(V^) 2 -m 2 0? 



d 2 v 



(11) 



Carrying out a Legendre transformation, it is easy to write down the Hamilto- 
nian density 



d 2 V 



where 



We also have 



2 V 



SL 3 a 
Pi = — = a A 4>i. 



2\ ( j )i (j) +\^ 2 k -v 2 )8 ij . 



Assume that the fluctuation field $^ decomposes into its constituents as: 



(12) 
(13) 
(14) 



=< $ > 



I* \ 

00 



(15) 
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The physical fields are defined so that 



Analogously, we define the classical background fields as: 



1 / 

"+ = ■^(vi+iV2); 



(16) 



(17) 



following the identification: 



( ^ \ 

V3 = v 

\ ° ) 



=< $ > 



(18) 



It is easy to see that a Legendre transformation provides a Hamiltonian given 



as: 



H — H-neutral ~t~ H charged "i" H-mixedi 



where 

Hneutral 



I 



(19) 



2a6 00 



+ 



(P 



1 



charged 



J d 3 xdta 3 (-^{(p 0+ po-) 



+ 1(V^)(V0 + ) + K)(0+0-) + ^(0^ ± - c4 ± )(0 + 0_)fr21) 



# m „ ed = J d 3 xa 3 dt{\(v 2 + 4> 2 _+v 2 _<p 2 + ) 
Where we have also put 



Po <Po 



A[2- 



,2l. ^ ^± 



A[2v+v_]; 



and 



2v + u_ + U3 + a 2 =v 2 . 



X[2a 2 } (23) 



(24) 



This is the most general Hamiltonian for the 0(4) scalar field system, now bro- 
ken up into a charged scalar field and two neutral scalar fields in the background 
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field formalism. The background field can now be parametrized through three 
angles: 

/ vCos{p)Sin{0)Sin{a) \ 
_ vCos(p)Sin{9)Cos(a) . . 

vSin{p)Sin{9) ' [2b > 

\ vCos(9) J 

Here we shall simplify the parametrization of the background field to two angles, 
9 and p by letting a = j: then, v± — -^Cos(p)Sin(9), V3 = vSin{p)Sin{9), 
and a = vCos(9). Canonical quantization gives the mode decomposed 
Hamiltonian: 



f d 3 k 1 ui^ t t 

Uneutral - J 

+^|(44 + 44) + _ i)(4 4 + 44 + d _ fcdfc + dt_ fcd t )} (26 ) 



H char ged = J ^3 (6 fc &/c + CfcC fc ) 



(2tt) 2 

~ 1)(&tfc5fc + cfectfe + 6 - fcCfe + ct -^ )} (27) 



/" d 3 fc Aa 3 w 2 cos 2 (p)sin 2 (0) t 

H mixed = / , ,o { : {bkb-k + bkcl 

J (2tt) 3 4o) 

+4&k + c fe c -fe + 4 c -fc + c fe fet fe + 6 I c fe + b l b -k) 
\a 3 v 2 cos{p)sin{p)sin 2 (9) 

(bka-k + ha{ + c\a k + c k a— k + c\a}_ k + c k a\ + b\a k + b\a}_^j 

\a 3 v 2 sin(p)sin(9)cos(9) / , + ,+ ,+ + \ 
— d k a- k + d k a\ + d\a k + d\a}_ k 

Xa 3 v 2 cos(p)sin(9)cos(9) 



+ 

+ 



(& fc d-k + 6 fc 4 + c{d k + c k d-k + c\S_ k + Cfe4 + &t d * + 6 I d -fc) > ( 28 ) 
where 

^l s <m = <m = {mi+ l ^ = (m|+ 4, (29) 

a 6 a 6 a 6 ^ a 2 a 6 a 2 
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and 

(I (It (JL (1) 

An important point to note here is that while ^l^(k) and oj^,(k) are momentum 
dependent, for ease of notation we will drop the k dependence for further calcu- 
lations and revive it when necessary. Clearly, there are two interesting cases to 
be considered here: 9 = and p — ^ . For the purposes of this study, we shall 
only consider the first case, 9 = 0. Its easy to see that H reduces to 

H = J -^ { 2^ { < ak + ak < ] 

~ 1 ^ a *° fc + afeatfe + a - kCLk + °-*4) + 2^(4 d * + d *4) 



4a 3 w| 
2a 3 V 2 



+ Cfc4) + - l)(b\b k + CfcC t + b _ kCk + c t_ fe6 t )} (31) 



This Hamiltonian has an sw(l, 1) symmetry and can be diagonalized by a 
scries of Bogolyubov transformations: In the neutral sector, 

A k (t, r) = ,i(r, t)a k + v(r, t)al k = U~\r, t)a k U{r, t). (32) 

Similarly for the sigma field, D k (t,r') and D\(t, r') arc defined in analogy with 
the definition of A k (t,r) and A k (t,r) with the d's replacing the a's. For the 
charged sector 

C k (t,r) =»c k + vbl k = U- 1 (r,t)c k U(r,t) (33) 
B k (t,r) = iiC- k + vb\ = U- 1 {r,t)b k U{r,t) (34) 

where 



A* = Cosh(r) = J1A + %.) + !] „ = Sinh(r) = J\[fi* + £*) - 1] 

and r is the squeezing parameter, satisfying ^ 2 -o 2 = 1 as required for a squeez- 
ing transformation. The complete unitary matrix for the squeezing transforma- 
tion may be written down as 

U(r t) = e J ll^ r ( fc ^){( a l at -A : - afca - fc )+(4 d -fc- d fc d -fc)+( c fc b -fc+ b fc c -fc)-( c I'' t -fc+''l c t J> 

(36) 
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Collecting all our results for the neutral and charged sectors, the total diag- 
onalizcd Hamiltonian is written in terms of various creation and annihilation 
operators as: 

/^3 u 1 1 1 

^^{^A(4^+^HCtC k +BtB k +l)}+n s (DlD k +-)}. (37) 

Since the S field decouples, we drop all terms associated with it whenever it is 
not essential to our arguments. The total dynamical Hamiltonian for the neutral 
and charged scalar fields in terms of the creation and annihilation operators 
(a, a\c,(J ,b and ^simplifies to: 

Ho - J^^{^ 2 + - 2 Mc k + blb k + l) 

+ iiv{{c k b- k + b k c_ k ) + {c{b ] _ k + b\cl k )} 

+ ^r{<y + v 2 ){a\a k + l) + vn{cJ_ k a\+a k a- k )}. (38) 

We have completed the derivation of a Hamiltonian describing the dynamics 
of neutral and charged scalar fields with a symmetry breaking parameter in 
an expanding FRW metric. We shall see that these two parameters allow us to 
examine the amplification of charged and neutral modes as they are varied, both 
in a sudden quench and adiabatically. We have already seen in [[5]] that this 
Hamiltonian can be used to describe the formation and decay of the disoriented 
chiral condensate (DCC). The generic structure of the Hamiltonian allows us 
to apply its dynamics to the inflationary studies in cosmology. We also believe 
it provides an excellent framework for analyses of the Bose-Einstcin condensate 
(BEC)- this is left for a later communication. 

3 Evolution of the fluctuations 

We note first that the sigma field decouples in this particular Hamiltonian and 
therefore it can be analyzed independently of the 4>o,4>± fields. We write the 
total dynamical Hamiltonian in terms of the creation and annihilation operators 
(a, a\ c, ,b and &^)and the squeezing parameters as: 

H = /(0^ 2 + ^ C ^ + 5 ^ + 1 > 
+ iiv{{c k b^ k + 6 fc c_ fe ) + (c\bl k + b\cl k )} 

+ -^{(^ 2 +^){a{a k + l} + (u f i){al k al + a k a- k }}. (39) 

Defining the following bilinear operators: 

V = a k a- k + b k c- k + c k b- k = Ky + + 
& = al k al + cl k b{ + bl k cl = K+ + K+ + K+ 
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N = \{ a \ a k + a^a-fe + b\b k + tf_ k b- k + c\c k + cl fe c_ fe + 3} 
= K° + K$ + K° 
it is easy to see that they satisfy an su(l, 1) algebra 

[AT, V] = -V; [N^] = D*; [2? f ,2?] = -2iV. (40) 
The su(l, 1) invariant Hamiltonian for the </>o, 0± fields assumes the form 

/d 3 k 1 
— — 2fy(fc, t)(^ 2 + ^ 2 )iV + 2fy(fc, i)/u/(X> + P f ) (41) 
(Z7r)° a° 

The time dependent evolution equation is given by 

H(t)\m>=ij t \m> (42) 
The particular su(l, 1) structure elucidated above provides us the solution: 

|^(t) >= e-^ (s£ r * (I> *- 1? * ) |^(0 > (43) 

for the evolution of the wave function immediately. Here r k is the squeezing 
parameter related to the physical variables fl^(k,t) and ui^(k) through 

, f*j,(fc,t) \2 _ j 

Where £l,p(k,t — ► oo) = Lu^(k). It is assumed that past and future infinity 
is flat spacetime. Thus in the evolution of the fluctuations, it is the frequency 
changes which bring about squeezing. 

The diagonalizcd Hamiltonian Hq can be converted into a Hamiltonian in 
terms of quantum fields corresponding to the operators A, B, C and their ad- 
joints to obtain a purely quadratic Hamiltonian. 

We can, for example, write 



n 



'j{A\A k + h = ^(A{A k + A k A{) = (% 2 U\(k,t) + P* A (k,t). 



Similarly for B and C [[5]]. The Hamiltonian H n can then be written as: 

d 3 k 

i=A,B,C 

The Schroedinger equation for each momentum mode is simply: 



(45) 



(46) 



H {k,t)^{k,t)=ij^{k,t). (47) 
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If we use the II-represcntation (co-ordinate space representation) for ip(k,t), 
then, the su(l,l) symmetry of the Hamiltonian tells us that the solution for 
ip(k, t) is just a Gaussian. The equation satisfied by the wave functions for each 
mode are then given by: 



$ A (M) + —4a + (^-) 2 (k,t)^ A (k,t) = 0. (48) 



3d ,• ,Q 
a "a 3 

(similar equations hold for fields B and C) where 

{ S hf { k,t) = (^)+m 2 +Xv 2 . (49) 

The expectation values of the number operator for the neutral scalar fields for 
each momentum k is given by: 

< Mt)\al a k\Mt) >= Smh 2 (r) =< ij> k \J$ k (t)A k (t)\il) k > ■ (50) 

An identical expression holds for the charged scalar fields. We let 

df] = a(i) _1 di (51) 

so that the FRW metric is transformed into a conformally flat metric: ds 2 = 
a(r)) 2 (dr] 2 — dx 2 ). The equations of motion given above are then transformed 
into ones that resemble a harmonic oscillator with time dependent frequencies. 
We shall write only the generic form of the above equations: In terms of the 
scaled time, rj, we have: 

In' 

i>" + + (k 2 + (m 2 + A(< $ 2 > -v 2 ))a 2 i> - (52) 

where a prime denotes differential wrt r\. 

Lastly, scaling if>: £ = aip so that the equation becomes: 

-i" + Vm = (k 2 +ml)£, (53) 

where 

V( V ) = a -!0 + m 2 (1 - a 2 ) - A(< $2 > (54) 
Thus in the symmetry broken stage < <f> 2 >= v 2 



V h U,) = a- 1 ^+ml(l-a i ) (55) 



and in the symmetry restored stage, < 3> 2 >= 

V r (r,)=a- 1 ^+ml(l-a 2 ) + \a 2 v 2 . (56) 

These equations can be interpreted in two ways allowing for simple solu- 
tions. Treating r\ as a spatial variable allows us to treat them as Schroedingcr 
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like equations with E — This then allows calculation of the reflection 

and transmission coefficients over the "potential barrier" provided by the V{rf) 
term. On the other hand, they are equations for time dependent harmonic oscil- 
lators with time dependent frequencies given by f2? and w?. These two pictures 
allow the calculation of the squeezing parameter dependent number operator 
N{k) = Sinh 2 (r(k)) involved in the evolution of the neutral and charged quan- 
tum fluctuations in the FRW background. 

We shall now look at the scenario in the FRW metric with expansion included 
to show that the enhancement of the low energy modes and the squeezing pa- 
rameter are dependent on the rate of the expansion mechanism by which the 
symmetry is restored. We have also seen that to produce substantial squeezing, 
we require a quenching. To show this really happens in this case, we have to 
compare the situation of a sudden quench with a slow adiabatic relaxation of 
the system from the symmetry restored stage to the symmetry broken stage. 

To examine the amplification of the neutral and charged modes for the case 
of a smoothly expanding metric, we shall assume the expansion parameter a{-q) 
to be[l]: 

a(r]) 2 = A + BTanh[brf\. (57) 

Here i measures the rate of the expansion and we choose b to be in units 
of m^. When viewed as a time dependent oscillator equation in an expanding 
metric we may write the equation for £ in the symmetry restored phase as an 
oscillator equation 

Z"+<4Z-V r {ri)Z = (58) 
and in the symmetry broken stage as 

?'+wl£-V b (r,)t = 0. (59) 

Where Vf,(?7)and V r (r]) are given by (56,55) with the particular choice of the 
expansion parameter a(r]) given by (57). Thus the change in frequency from the 
restored to the broken stage is given by 

V r (r ) )-V b (r ) ) = +\v 2 a(r 1 ) 2 (60) 

Here, we shall choose A and B in eqn.(57) to be 1. The transmission coeffi- 
cients across such a barrier are easily looked up from [1]. While the number of 
particles of mode k is given by: 



{Sinh 2 {^-)) 



Here b measures the duration of the quench and its effect is exhibited in figure 
1. 

N(k)vs.k is plotted in fig. 2. We see that in the adiabatic limit of large 
6, N(k) is exponentially suppressed so that there is no enhancement of low 
momentum modes. 
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Figure 1 : Shows the transition from quench to adiabatic for various values of b 
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Figure 2: Shows the variation N(k) with k for values of b=0.35 (red); 0.4 (green) 
0.5 (blue). 
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From the above we conclude that if expansion is rapid we get the quenched 
limit, while, if expansion is slow, we get the adiabatic limit. Since in the squeezed 
state description N(k) = Sinh 2 rk , where = r(k) is the squeezing parameter , 
we conclude that in the quenched limit the squeezing parameter is much greater 
than in the adiabatic limit. This demonstrates clearly the connection between 
the rate of expansion and squeezing and the enhancement of charged and neutral 
particle creation. 

Since the amplification is large for the low momentum modes, let us look at 
the charged and neutral number particle distributions for zero momentum. We 
shall add a subscript to denote zero momentum modes so that the distributions 
are given by: 

Pn ,n + ,n- = | < no, n+, rc_ |V> > | 2 

= < n |e r °( a ^ 2 - ro * a « |0 >< n+, rc_ | e 2r °( 6 H-*>oco)| > |2 ^) 

defining S(ro) as the one mode squeezing operator 

S(ro) =< n |e ro «^ )2 - a o)|0 >= S no , . (63) 
S tm (ro) is then the two mode squeezing operator 

< n + ,n_|e (ro(b « c ^ &oco)) |0 >= S^ n _ a . (64) 
The neutral and charged particle number distribution is: 

Pn„,n c —< S no fi > 2 < S^™ + n _ a o > 2 (65) 

Writing n + = n_ = n c , we get the distribution of charged particles to be 

^ (tonfe(ro)) 2 "- 
P nc =^P no , nc = {cosh{ro))2 (66) 



and the distribution of neutral particles to be 

P _Y- P n \(tanh(r )) n ° 



Thus the neutral and charged particle number distributions are significantly 
different as the two types of squeezed states appearing in the expressions above 
have different properties. We now illustrate the effect of squeezing in these two 
distributions. Figs. 3, and 4 show the difference in the charged and neutral 
particle distributions as we vary the squeezing parameter from a low value to a 
high value. 

We now illustrate the effect of quenching versus adiabaticity on these two 
distributions. Figs. 5 and 6 show the difference in the charged and neutral 
particle number distributions as we vary from the adiabatic limit where the 
difference is negligible to the quenched limit where the difference is significant. 
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Figure 3: Shows the variation of P nc with n for the r = 3, 3.5 and 4 
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Figure 4: Shows the variation of P nc with n for the r$ = 3, 3.5 and 4 




Figure 5: Shows the variation of P„o (solid line) and P nc (dashed line) with n 
the adiabatic limit (r = 2) 
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Figure 6: Shows the variation of P„o (solid line)and P nc (dashed line) with n for 
the quenched limit (r = 4) 

4 Conclusion. 

To conclude, in this paper we have constructed an effective Hamiltonian for the 
study of the dynamics of charged and neutral scalar particles in an expand- 
ing background FRW metric. Starting from an 0(4) sigma model through the 
inclusion of second order quantum fluctuations we have shown that the Hamil- 
tonian is quite generic and in addition to applications to the formation of the 
DCC [5] , can also be used to study the enhancement of low momentum charged 
and neutral particle modes in an FRW metric suggesting an application to re- 
cent ideas of "preheating" in an inflationary cosmological scenario. We have 
seen the appearance of SU(1,1) symmetries leading to the presence of squeezed 
states in their dynamics. We find that in the quenched limit ( fast expansion) 
the low momentum modes are enhanced significantly, whereas in the adiabatic 
(slow expansion) limit, no such enhancement occurs. This has been shown to 
be directly related to the value of the squeezing parameter. The manifestation 
of this difference shows up directly in the total neutral and charged number 
distributions at zero momentum. 
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